Landau theory of charge and spin ordering in the nickelates 
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Guided by experiment and band structure, we introduce and study a phenomenological Landau 
theory for the unusual charge and spin ordering associated with the Mott transition in the perovskite 
nickelates, with chemical formula RNiOs, where R=Pr, Nd,Sm, Eu, Ho, Y, and Lu. While the 
Landau theory has general applicability, we show that for the most conducting materials, R=Pr, Nd, 
both types of order can be understood in terms of a nearly-nested spin density wave. Furthermore, 
we argue that in this regime, the charge ordering is reliant upon the orthorhombic symmetry of the 
sample, and therefore proportional to the magnitude of the orthorhombic distortion. The first order 
nature of the phase transitions is also explained. We briefly show by example how the theory is 
readily adapted to modified geometries such as nickelate films. 



The nature of the Mott metal-insulator transition 
(MIT), driven by Coulomb repulsion between electrons, 
is a central subject in condensed matter physics. A 
canonical MIT occurs in the nickelates, perovskites with 
the composition RNiOa, where R is a rare earth with 
nominal valence R^+. Metallicity in the nickelates cor- 
relates with ionic radius, varying from largest, R=La, 
which is metallic at all temperatures, to smallest, R=Lu, 
which is insulating at all temperatures • R=Eu has the 
highest MIT temperature, Tmit 480K.fflThe nickelates 
are particularly interesting because they display complex 
ordering phenomena in the insulating state. The insulat- 
ing ground states are magnetic, with a large unit cell cor- 
responding to a periodicity of 4 lattice spacings relative 
to the ideal cubic structure, and is usually interpreted in 
terms of an "up-up-down-down" spin configuration. Re- 
markably, this complex ordering seems to be consistent 
across the entire family.l^H^ This magnetism coexists with 
a form of "Charge Ordering" (CO), which modulates the 
Ni charge from Ni^^"^ to Ni^""^ on alternating sites of a 
rock-salt type substructure of the cubic perovskite. In 
the more insulating nickelates, R=Eu, Ho, this charge 
ordering occurs not only in the ground state but also 
in an intermediate temperature insulating phase without 
magnetism between the high-temperature metallic phase 
and the low temperature magnetic one. 

It has been postulated that the charge ordering is fun- 
damental, and should be thought of as separation into 
spinless Ni^+ and spin 5 = 1 Ni^+ ions, due to dom- 
inant local Hund's rule coupling Jh^ However, no ex- 
planation has been given for the particular complex yet 
robust magnetic structure in this picture. Moreover, in 
the materials PrNiOs and NdNiOs, for which the MIT 
occurs at relative low temperature, the ionic description 
does not seem so natural. In addition, in these materials, 
charge and spin order occur simultaneously, bringing the 
primacy of the former into question. 

In this paper, we introduce a phenomenological Lan- 
dau theory, motivated by the particular structure of the 
e^ bands of the nickelates near the Fermi energy. It is 
particularly relevant to the more itinerant materials such 
as PrNiOs and NdNiOs, where charge and spin ordering 



are coincident, and the Landau theory has a natural mi- 
croscopic origin in a spin density wave (SDW) picture. 
The SDW picture and Landau theory yield a number of 
direct insights into the ordering in bulk samples. First, 
the observed unusual magnetic periodicity is simply ex- 
plained due to approximate Fermi surface nesting of the 
Cg bands. Second, we argue that the observed charge or- 
dering is naturally induced in this picture as a secondary 
order parameter. Third, the magnitude of the induced 
charge order is proportional to the degree of orthorhom- 
bicity in the material, a correlation which is indeed ob- 
served in experiments, though the causative relationship 
does not appear to have been identified previously. 

The Landau theory can also be easily adapted to other 
situations. As one such application, we derive the a gen- 
eral phase diagram describing the passage from "weak" 
(SDW-like) Mott insulators to "strong" Mott insulators, 
in which charge ordering is indeed dominant, which is 
in agreement with that observed in the bulk nickelates, 
and we explain the first order nature of the transitions 
observed there. We also show how to incorporate the ef- 
fects of strain and interfaces in films and multilayers, and 
predict that changes of the magnetic ordering wavevector 
and exotic "multiple-q" states may be thus induced. 

Microscopic considerations: It is helpful to consider 
as a semi-microscopic model framework a tight-binding 
description of the eg bands, 

Htb = -Y.tt!cL(^j,., (1) 

where i,j are site indices, a, 6 = 1,2 are orbital indices 
denoting the two eg states with the symmetry of — 
x'^ — y'^ and x'^ — y'^ orbitals (these should properly be 
considered hybridized combinations - Wannier functions 
- of Ni d states and the neighboring O p states, as RNiOs 
is a charge transfer insulator), and a =t, i is the spin 
index. Sums over orbital and spin indices are implied. 
The dominant hopping processes are expected to be those 
with cr-type bonding, which can occur with amplitude 
t and when z,j are first and second neighbor sites, 
respectively, of the ideal cubic Ni sublattice. Specifically, 
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FIG. 1: Fermi surfaces for the tight-binding model. In (a) and 

(b) , we show the conduction and valence band Fermi surfaces, 
respectively, for t' /t = 0.1. For larger t' /t^ the conduction 
band Fermi surfaces become large and hole-like, as shown in 

(c) and (d) for t' /t — 0.3. The approximate nesting in the 
latter case is indicated schematically in (d). 
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FIG. 2: Zero frequency spin susceptibility for the tight- 
binding Hamiltonian for t' /t — 0.1,0.2,0.3, as a function of 
momentum k in the cubic Brillouin zone. Note that for the 
best nested situation, t' /t — 0.3, the susceptibility is sharply 
peaked close to the wavevector 27r(|, |, \). 



then t^^^ = ^M±A±^ = ^'(^m^' + ^mO' 

where = = ("i-^), = (1,0) are 

the orbital wavefunctions for the 2x^ —y'^ — z^^ 2y'^ — x^ — 
z^, and 2z^ — — y'^ a-bonding orbitals along the three 
axes, in the chosen basis, respectively. 

The above tight-binding model agrees with the eg 
bands obtained from LDA calculations,^ and the Fermi 
surface measured recently in photoemission on a LaNiOa 
film.^ It also explains well resistivity. Hall effect, and 
thermopower measurements on LaNiOa films, ^ and ap- 
proximately describes the inter-band optical spectral 
weight at energies below about leV.^^ LDA favors t' /t ^ 
0.1, while comparison to photoemission is best for t' /t ^ 
0.3, and the range 0.1 < t' /t < 0.4 is consistent with 
transport. It is instructive to view the Fermi surface for 
this range of values (see Figjl]). One observes in the mid- 
dle of this range that the Fermi surface is rather flat and 
approximately composed of "cubes" rather than spheres. 
Fermi surfaces with large flat sections are considered ap- 
proximately "nested" , and well-known to lead to enhance 
susceptibilities at certain wavevectors. Calculation of the 
free electron spin susceptibility (see Fig. [2| indeed shows 
an enhancement, peaked about wavevectors (kkk)^ with 
k close to |(x27r) (for t^t = 0.3 we find e.g. k ^ OAir). 
With interactions included, the random phase approx- 
imation, or Hartree-Fock theory, both thereby show a 
tendency to spin- density- wave (SDW) order at wavevec- 
tors close to this one. It is remarkable that band theory 
considerations give a simple mechanism for SDW order at 
k = 7r/2, which corresponds precisely to that observed in 
all the insulating nickelates. Interestingly, we have also 
found a mechanism for this order in the strong-coupling 
limit of the appropriate multi-orbital Hubbard model, 
which is too involved to report here.^^ Together, this may 
explain the robustness of the magnetic state observed in 
experiment. This unusual magnetic order has not, to our 
knowledge, been explained before by any theory. 

Landau theory: Rather than proceeding with a micro- 
scopic theory (Hartree-Fock and other calculations will 
be reported in a future publication^) , we instead pursue 
the implications of this view using symmetry-based Lan- 
dau analysis. We begin by considering the problem for 
an ideal cubic solid, and take into account distortions of 



the perovskite structure at a later stage. For the SDW 
order, states with wavevectors along any of the (111) 
axes are equivalent. Hence we actually need to include 
four order parameters, with wavevectors Qa^ given by 
Qo = 1(1,1,1), Qi = 1(1,-1, -1),Q2 = f(-l,l,-l), 
and Qs = f(— 1,-1,1). Physically, the meaning of the 
order parameters, which are complex vectors, is that 

(S,) = ^Re [Vae*^-'--] . (2) 

a 

It is a straightforward but lengthy exercise to identify 
all the allowed terms, up to fourth order in the SDW 
order parameters, in the Landau expansion, taking into 
account the translational and point group symmetries of 
the ideal cubic structure. At quadratic order, one obtains 
only the single coefficient 

F, = rJ2ra-^a^rY,\^a\^- (3) 

a a 

At quartic order, one finds 13 distinct terms. We divide 
them into the three terms involving products of only one 
"flavor", 

a a 

+U3(^(V„-V'„)' + h.c.), (4) 

a 

and ten additional terms involving products of two and 
four distinct fields. For brevity, we will not give these 
explicitly here, as they will not play a major role in what 
follows. 

Single-q states: The full free energy simplifies greatly 
if we restrict to "single-q" states, consistent with experi- 
ment. Here only one of the four order parameters is non- 
zero, and the ten quartic Landau invariants not shown 
in Eq. Q vanish. In a single-q state, we denote below 
the non-zero order parameter by simply '0, without a 
subscript. 

Several types of single-q states are possible, dependent 
upon values of U2 and us in ^4^^ in Eq. Q. The second 
coefficient, 2/2, distinguishes between collinear and spiral 
states. For U2 > 2\us\^ minimum energy states have ijj = 
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'0(ni+m2), where lii and n2 are orthogonal unit vectors. 
These describe coplanar spirals, with spins of fixed length 
in the plane spanned by lii and n2. The phase of is 
arbitrary, corresponding to rotations of the spins within 
this plane. 

For U2 < 2\us\, the free energy is minimized by i/? = 
V^fi, where -0 is a complex scalar. Such configurations 
describe, via Eq. ([2]), collinear SDW states. In this case, 
the remaining coefficient, u^^ selects preferential phases 
of ip. Writing ip = IV^Ie*^, we finally have 

Fi'^ = K + U2M^ + 2us\^\^ COS40. (5) 

When < 0, states with = are favored, while 
for Us > 0, states with = ^{n -\- ^) are favored (in 
both cases n = 0, 1, 2, 3 describe four spatially translated 
states). Going back to Eq. ([2|, and using the wavevector 
Qq = 1^(1,1,1) and taking n = z for concreteness, we 
find that these describe spin states with 

1^1 X (+1,0,-1,0,---) forw3<0 

(6) 

where the successive terms in parenthesis describe succes- 
sive spin expectation values when moving in unit steps 
along the principle cubic x, or z axes in real space. 

Charge order: For the three types of magnetic states 
found above (one spiral, two collinear), let us consider the 
associated charge order. We consider charge ordering at 
the wavevector (tt, tt, tt), corresponding to the "rock salt" 
ordering observed in experiment. The order parameter 
^ is introduced via 

(n,)=n + (-l)^^+^^+^^^, (7) 

where = ^laa^iaa electron number oper- 

ator. $ may equally well be regarded as representing 
the amplitude of optical phonons representing octahedral 
breathing modes at the same wavevector. Symmetry al- 
lows the following terms in the free energy involving ^ 

F$ = + ^$4 _ A$ ^ Re [i,^ ■ -0,] , (8) 

a 

where we have included the leading linear coupling to the 
SDW order parameters. The crucial thing to note here is 
that whenever i/^ • i/? is real and non-zero, a non-zero ^ is 
necessarily induced in the minimum free energy state. In 
such situations, $ is a secondary order parameter, slaved 
to the primary SDW one. Mathematically, if u can be 
neglected, ^ can be readily "integrated out" {F^ can be 
minimized with respect to $) to simply renormalize the 
quart ic SDW couplings. 

In this case, we can analyze charge order simply in 
terms of the SDW states. In two of the three SDW states 
discussed above, the spiral and the collinear state with 
6> = 7r/4, Re['0 • -0] =0, the charge order vanishes. This 
is easily understood since in these cases all sites have 
equivalent spin states, up to rotations. In the remain- 
ing collinear state, with = 0, Re['0 • -0] = and 



^ 7^ 0. Again this is intuitively clear since the sites with 
zero spin are obviously distinct from those with non-zero 
spin. Two problems arise now in comparison with ex- 
periment. First, we must arbitrarily choose parameters 
to be in the third magnetic state in order to obtain the 
observed charge order. Second, in experiment, non-zero 
magnetic moments, of unequal magnitude, are clearly ob- 
served on the two types of "inequivalent" sites, while in 
this theoretical configuration the magnetic moment van- 
ishes on half the sites. 

Orthorhombicity: These problems are resolved by tak- 
ing into account the distortions of the ideal perovskite 
structure. We focus on the orthorhombic (Pbnm) struc- 
ture, which obtains for all the nickelates except LaNiOa, 
which is rhombohedral, and does not undergo a Mott 
transition. The orthorhombic distortion is present in 
the metallic state up to high temperatures (e.g. up to 
T = 780K in PrNiOa^^), and is understood to arise from 
reduction in the tolerance factor due to the changing 
rare earth ionic radius. The Pbnm space group has only 
discrete reflection and inversion operations in its point 
group. The structure has a quadrupled unit cell, which 
is doubled by a 45° rotated, approximately ^/2 x en- 
larged supercell in the a — b plane, and a doubling along 
the c axis. One can rewrite the conventional cubic coor- 
dinates x^y^z'm. terms of standard orthorhombic coordi- 
nates X, y, z, according tox = x + y— ^,7/ = — x + y + ^, 
2: = 2z. Making this transformation, one finds that the 
4 cubic SDW states corresponds to two orthorhombic 
SDW wavevectors: Qg^*^° = Q§^*^° = 27r(0, ^, |) and 
gortho ^ gortho ^ 27r(^, 0, \) . The latter is the wavevec- 
tor found in experiments on nickelates with R=Sm, Eu, 
Nd, Pr, and Ho.^ ^ 

The Pbnm space group has considerably lower sym- 
metry than cubic, containing only inversion, reflection, 
and 180° screw axes apart from translations. As a con- 
sequence, the orthorhombic distortion allows additional 
terms in the Landau free energy. A straightforward anal- 
ysis shows that two such terms arise at quadratic order: 

= ri(|^op-|^ip-|^2p + |^3p) 

+r2 (V^i • V^i - '02 • '02 + c.c.) . (9) 

We neglect orthorhombic corrections to the quartic terms 
on the grounds that they are presumably smaller. 

The terms in Eq. ([9| clearly lead to an energy differ- 
ence between the two possible orthorhombic wavevectors. 
We focus on the case ri > — |r2|, for which ordering in 
01, 02 is preferred, which corresponds to the experimen- 
tal wavevector. To proceed, we suppose, without loss of 
generality, that 0i = is the non-zero order parameter. 
For the spiral state, we now consider the cubic energy 
from Eq. Q, with ii2 > 0, added to the r2 term. For r2 
non-zero, the minimum energy configuration is deformed 
to (X (1 + 5)fvi + i{l — (5)n2, with (5 oc r2. This corre- 
sponds to a deformed spiral in which spins trace out an 
ellipse rather than a circle. In this state, <I> oc • ex r2, 
so charge order is induced. 
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between collinear and spiral spin states, taking i/? = ^^x 
with ip real. Then the full Landau free energy is 



FIG. 3: Mean- field phase diagram for the Landau theory of 



Eq. (11), describing the interplay of spin and charge order- 
ing. The dashed (solid) lines delineate first (second) order 
transitions within mean field theory. 



For the collinear state, Eq. ([5| is modified to 

F{0) = Fo + 2^X3|V^|'^cos4(9 + 2r2|V^pcos2i9. (10) 

There are two cases to consider here. For < 0, both 
cosines can be simultaneously minimized. The minimum 
in this case occurs at either = ^n^ with integer n either 
even or odd if r2 is negative or positive, respectively. This 
selects the +1, 0, — 1, 0, • • • type ordering in the first line 
of Eq. ([6|. 

In the other case, 2x3 > 0, the situation is more inter- 
esting. The minimum at ^ = |^(n + 1/2) are unstable to 
small r2 of either sign. Thus the generic situation in this 
case, at least for small |r2|, is to obtain a generic value of 
< ^ < 7r/4. (For sufficiently large r2, the minimum 
will again lock to the above type of solution). As a result, 
one obtains the type of order observed in experiment. 

The conclusion is that in the = 27r(^,0, ^) 

state, the SDW is always accompanied by charge order 
due to orthorhombicity. Moreover, if the state is such 
that moments are present on all sites, then the charge 
order is proportional to the degree of orthorhombicity, 
parametrized here by r2. This is an important conclusion 
of this paper. Interestingly, if the analysis is repeated for 
a rhombohedral (R3c) crystal, no charge order is induced 
in this way.^^ Hence for a rhombohedral insulator, we 
predict the occurrence of a "pure" SDW state. 

Order of transitions: The three phase transitions - 
metal-CO, metal-SDW, and CO-SDW - observed in ex- 
periment are all first order. This behavior can be ratio- 
nalized by various mechanisms. First consider a mean- 
field treatment of the Landau free energy. For simplicity, 
we assume single-q states, and neglect the distinctions 



F = rilj^ ^ ui/j^ - X^i/j^ + r^^ + u^^. 



(11) 



By minimizing this free energy, one obtains the phase 
diagram shown in Fig. [3] For A = 0, the SDW and 
CO orders are tuned independently by r and f respec- 
tively. However, with non-zero spin-charge coupling A, 
the region near the origin is modified. In particular, 
the consequence is that for systems close to the multi- 
critical point at which the CO phase emerges, all the 
transitions to the SDW phase become first order. This 
may explain the discontinuous SDW transitions seen in 
experiment. To explain the first order metal-CO transi- 
tion, we note that electron-phonon interactions are signif- 
icant in nickelates.^^ In particular coupling of ^ to acous- 
tic phonons renders the theory of the metal-CO transi- 
tion equivalent to the compressible Ising models which is 
known to have a first order transition.^ 

Films: We conclude with examples of how the Landau 
theory is adapted to modified geometries. Pure strain is 
relevant to thick films in which interfaces are not impor- 
tant. For a tetragonal substrate, symmetry implies that 
the strain modifies r (Eq. ([s])) and ri (Eq. (|9)). If the 
intrinsic orthorhombic contribution to ri is small, we pre- 
dict that strain can switch the wavevector from Q°y2^° to 

(5q^3^°, and at the same time "turn off" the charge order- 
ing. A more severe effect occurs at an interface, due to 
the lack of translational symmetry normal it. This leads 
to symmetry-allowed terms which mix the SDW states 
with wavevectors Qi and Q2 (or Qq and Qs)^ so that 
"multiple-q" ordering appears near the interface. 



Acknowledgments 

We are grateful to Susanne Stemmer, Jim Allen, Dan 
Ouellette, and Junwoo Son for discussions and exper- 
imental inspiration. This work was supported by the 
NSF through grants PHY05-51164 and DMR-0804564, 
and the Army Research Office through MURI grant No. 
W911-NF-09-1-0398. 



^ J. Torrance et al., Phys. Rev. B 45, 8209 (1992). 

^ M. T. Fernandez-Diaz et al., Phys. Rev. B 64, 144417 
(2001). 

^ J. Rodrfguez-Carvajal et al., Phys. Rev. B 57, 456 (1998). 

^ J. L. Garcia-Muhoz et al., Phys. Rev. B 50, 978 (1994). 

^ J. Garcia-Mufioz et al., Europhys. Lett 20, 241 (1992). 

^ I. Mazin et al., Phys. Rev. Lett 98, 176406 (2007). 

^ N. Hamada, J. Phys. Chem. Solids 54, 1157 (1993). 

^ R. Eguchi et al., Phys. Rev. B 79, 115122 (2009). 



^ J. Son et al., Appl. Phys. Lett 96, 062114 (pages 3) (2010). 
D. G. Ouellette et al., Phys. Rev. B (2010), submitted. 
S. B. Lee et al., in preparation. 

P. Lacorre et al., J. Solid. State. Chem 91, 225 (1991). 
^2 J. Vassiliou et al., J. Solid. State. Chem 81, 208 (1989). 
M. Medarde et al., Phys. Rev. Lett. 80, 2397 (1998). 
D. J. Bergman et al., Phys. Rev. B 13, 2145 (1976). 



